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ABSTRACT 



Transverse oscillations of thin threads in solar prominences are frequently reported in high-resolution observations. Two typical 
features of the observations are that the oscillations are damped in time and that simultaneous mass flows along the threads are 
detected. Flows cause the dense threads to move along the prominence magnetic structure while the threads are oscillating. The 
oscillations have been interpreted in terms of standing magnetohydrodynamic (MHD) kink waves of the magnetic flux tubes which 
support the threads. The damping is most likely due to resonant absorption caused by plasma inhomogeneity. The technique of 
seismology uses the observations combined with MHD wave theory to estimate prominence physical parameters. This paper presents 
a theoretical study of the joint elfect of flow and resonant absorption on the amplitude of standing kink waves in prominence threads. 
We find that flow and resonant absorption can either be competing effects on the amplitude or both can contribute to damp the 
oscillations depending on the instantaneous position of the thread within the prominence magnetic structure. The amplitude profile 
deviates from the classic exponential profile of resonantly damped kink waves in static flux tubes. Flow also introduces a progressive 
shift of the oscillation period compared to the static case, although this elfect is in general of minor importance. We test the robustness 
of seismological estimates by using synthetic data aiming to mimic real observations. The effect of the thread flow can significantly 
affect the estimation of the transverse inhomogeneity length scale. The presence of random background noise adds uncertainty to this 
estimation. Caution needs to be paid to the seismological estimates that do not take the influence of flow into account. 

Key words. Sun: filaments, prominences — Sun: oscillations — Sun: corona — Magnetohydrodynamic s (MHD) — Waves 



1. Introduction 



filaments are large-scale magnetic 
corona ( see th e recent reviews by 



Solar prominences and 

structures of the solar 

iLabrosse et al .'2010; Mackav et a l. 2010, about the physics, dy 
namics, and modelling of prominences). High-resolution Ha 
observations reveal that prominences are formed by a myr- 
iad o f fine structures usually called threads (see, e.g., iLinI 
1201 Ih . Threads are thin and long plasma condensations which 
outline the prominence magnetic field. Theoretically, promi- 
nence threads have been modelled as magnetic flux tubes an- 
chored in the sola r photosphere (e.g.. iBallester & Priest! [19891 
iRempel et al.l l999). which are only partially filled with the cool 
(~ 10"* K) prominence material, i.e., the thread itself, while the 
rest of the tube is occupied by hot coronal plasma and therefore 
invisible in Ha images. 

Threads are highly dynamic. For example, transverse thread 
oscillations and propagating waves along the threads are fre- 
quently observed in prominences, which have been interpreted 
in terms o f magnetohydrodynamic (MHD) kink waves (se e the 
reviews bv lArregui & Ballesteiir2010l:lArregui et al.ll2012ah . The 
reported periods are usually in between 2 and 10 minutes, 
and the osc illations are typically damped after a few periods 
dNing et alj 12009.) . Resonant absorption, caused by plasma in- 
homogeneity in the transverse direction to the magn etic field has 
been proposed as the damping mechanism (e.g., Arregui et alj 
[2008, 2011; Soler etal. 2009, 2010.) . In addition, flows and 
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mas s motions in prominences have been also rep orted (e.g., 
ZirkeTetaD \l99^ iLin et all [200l lOkamoto et aTT r2007). The 
mean flow velocities are in most cases less than 30 km s"', al- 
though values up to 40-50 km s"' have been observed in ac- 
tive region prominences. The presence of flow can have a direct 
impact on the behavior of the waves (see a discussion on this 
issue in[Carbonell et al. 2009). The combined use of both obser- 
vations and MHD wave theory allow s to apply the technique of 
seismology to solar prominences (see iOliveij2009t[Arregui et al.l 
[20T2bl) . 

Some observations point out the simultaneous presence of 
transverse oscillations and mass flows in prominence threads. 
The work by [O kamoto et al. ( 2007) is probably the best exam- 
ple. [Ok^otoetaL (2007) observed an active region prominence 
with Hinode/SOT using Ca II H-line images. They detected that 
some threads in the prominence were flowing presumably along 
the magnetic field with an apparent velocity on the plane of sky 
of around 40 km s"'. Simultaneously, the threads were oscillat- 
ing in the transverse direction. The mean period of the oscilla- 
tions was 3 min. The oscillations were in phase along the whole 
length of the threads, which was roughly between 3,000 km 
and 16,000 km. Thus, a standing wave interpretation was pro- 
posed and the actual wavelength was estimated to be at least 
250,000 km, which would correspond to twice the total length 
of magnetic field lines, approximately. 

The observations bv lOkamoto et al.l (l2007l) motivated a num- 
ber of the oretical works that aim to interpret the observed thread 
dvnamics. lTerradas et al.l (l2008l) interpreted the observations by 



1 



R. Soler et al.: Damped oscillations of flowing threads 



lOkamoto et alJ (l2007h in terms of standing MHD kink modes 
and used the observed periods to perform a seismological es- 
timation of a lower bound of the pro minence Alfven speed. 
Subsequentlv. lSoler & GoossensI (|201 ih revisited the same event 
and performed both an analytical and a numerical investiga- 
tion of the influence of flow on the period and amplitud e of 
the standing waves. However, neither iTerradas et al.l (1200 8*) nor 
[Soler & Goossens (201 1J took the damping of the oscillations 
into account. Okamoto et al.l (l2007h did not study the evolution 
in time of the amplitude of the transverse thread oscil lations, but 
damping has been reported in other events (see, e.g.. lNing et alj 
12009 ). The main purpose of the present work is therefore to de- 
termine the joint effect of flow and resonant absorption on the 
amplitude of standing kink MHD modes in prominence threads. 
The resonant damping of kink waves in the presence of flow has 
been studied in the past in the context of coronal loops (see, e.g.. 



iGoossens et al .ll992l:lTerradas et al.l2010l:ISoler et al.l201 ih . but 

not in the case of prominences. 

We apply th e general analytic theory developed by 
iRudermanlfeOllai ^ for standing kink waves in magnetic flux 
tubes with a time-dependent background. Analytic expressions 
for the kink mode amplitude and period as functions of the 
model parameters are derived. In addition, since the efficiency 
of resonant damping is directly controlled by the transverse in- 
homogeneity length scale of the threads, part of the present paper 
is also devoted to test the impact of flow on the e stimation of this 
para meter using the technique of seismology ( iGoossens et alj 
I2OO8 ). 

This paper is organized as follows. Section |2] contains a de- 
scription of the prominence thread model used in this work. 
The mathematical method is presented in Section [3] where the 
main equations are obtained. Then, the impact of the different 
model parameters on the amplitude of kink modes is studied in 
Section m Later, the implications of our results for prominence 
seismology are discussed in Section |5] Finally, the conclusions 
of this work are given in Section |6] 

2. Model 

The prominence thread model used in this paper is schematically 
shown in Figure [T] It is composed of a straight and cylindrical 
magnetic flux tube of radius R and length L. The z-axis is cho- 
sen so that it coincides with the axis of the tube. The magnetic 
field is B = Be,, with B constant everywhere. The tube ends 
are located at z = +L/2 and are fixed by two rigid walls rep- 
resenting the solar photosphere. The magnetic tube is partially 
filled with prominence plasma, i.e., the prominence thread, with 
length Lp and density pp, while the rest of the tube, i.e., the evac- 
uated part, is occupied by less dense plasma of density pe- The 
external plasma represents the coronal medium with density pc. 
According to the observed typical values of thread widths and 
lengths from the high-resolution observations of prominences 
(see, e.g., Lin 201 1), the ranges of realistic values of R and Lp are 
50 km < < 300 km and 3,000 km < Lp < 28,000 km. The total 
tube length, L, is very difficult to measure from the observations, 
but one can relate L to the typical spatial scale in prominences 
and filaments, i.e., L ~ 10^ km. These estimations of L, Lp, and 
R imply that prominence threads are very thin and long struc- 
tures. The ratio of the prominence density to the coronal density 
is a large parameter. The value Pp/pc - 200 is usually consid- 
ered in the literature, although we do not fix the density contrast 
at the present stage. For simplicity we assume that the evacuated 
part of the tube has the same density as the corona so that we 
takepe - Pc- 



In addition, in the prominence thread we consider a trans- 
versely inhomogeneous transitional layer of thickness I, that 
continuously connects the internal prominence density, pp, to the 
external coronal density, pc. The limits l/R = Q and IjR -2 cor- 
respond to a thread without transverse transitional layer and a 
fully inhomogeneous thread in the radial direction, respectively. 
The present paper adopts the /3 - approximation, where /3 
refers to the ratio of the gas pressure to the magnetic pressure. 
Hence, the plasma temperature is irrelevant for this study and we 
can choose the density profile in the transitional layer arbitrar- 
ily. The presence of this transitional layer introd uces damping of 
the kink waves by resonant absorption (see, e.g. . IGoossens et aTl 
I1992L I2OO6I I2OI II) . Finally, as a key ingredient of the model, 
there is a flow along the tube with a constant velocity uq. As a 
consequence, the prominence thread moves along the tube as a 
block at the velocity mq. 

Standing kink waves supported by the present model have 
been investigated in the past under several simplifications. In 
the absence of flow and for l/R — 0, ther e are the studies by 
iJoarder et a and 'Diaz et akl (l200ll) in Cartesian geom- 

etry, and by[Diazetal. (2002) in cylindrical geometry. These 
works considered arbitrary values of R and L, which involves 
a com plicated mathematical formalism. iDvmova & RudermanI 
(l2005l) derived a more simplified approach valid for thin cylin- 
drical tubes, i.e., RjL «: 1. The thin tube approximation 
provides very accurate results when compared to the gen- 
eral solution. Subsequently, the case IjR + in t he ab sence 
of flow w as addressed by IDvmova & RudermanI (12006 *) and 



Soler et al.l (120 10) again in the thin tube approximation, and by 
Arregui et aljl201 11) in the general case. In the presence of flow, 
the pr oblem was studied by^Terradas et al. (2008) numerically 
and bv lSoler & GoossensI ([201 lb combining analytic and numer- 
ical methods. However, both works considered IjR - and so 
resonant damping was absent from their investigations. 

3. Mathematical method 

The mathematical method followed in thi s work is based on 
the theory developed by RudermanI (1201 lalfbh . who derived the 
governing equation describing standing kink waves in the thin 
tu be and thin boundary approximations (see the Equation (23) 
of lRuderm an 2011b). We refer the reader to the original paper 
by iRudermanl (1201 Ibl) for a detailed derivation of the govern- 
ing equation. The thin tube approximation means that we are 
restricted to the case RjL <si 1 and R/Lp <sc 1. This is not a prob- 
lem since the realistic values of L, Lp, and R satisfy these con- 
straints. The thin boundary approximation means that we take 
l/R <s: 1, i.e., an abrupt but continuous transition in density in 

the transverse direction^ 

Equation (23) of IRudermanl ( 1201 Ibh can be solved analyti- 
cally using the Wentzel-Kramers-Brillouin (WKB) approxima- 
tion (see, e.g.. Bender & Orszag 1978). The WKB approxima- 
tion is appropriate when the background properties are slowly 
varying functions of space and/or time. In the present applica- 
tion of the WKB approximation we assume that the time scale 
related to the waves, e.g., the period, is much shorter than the 
time scale related to the changes of the background configura- 
tion. As in lSoler & Goosseni (1201 ih we define the parameter 6 



as 



6 = 



Mo 



(1) 



with Mo the flow velocity and L the total length of the tube. The 
validity of the WKB approximation is restricted to small values 
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of 6 so that P6 ^ I, where P is the period of the sta nding kink 
wave. In the observations by lOkamoto et alJ (l2007i) . the mean 
flow velocity and period are uq » 40 km s"' and P ^ 3 min. 
For L ~ 10^ km these values give Pd x 0.072. Hence the condi- 
tion of applicability of the WKB approximati on is fulfilled in the 
present case. Using the WKB approximation. lRudermanl(l201 Ibh 
arrives at the expression for the transverse displacement of the 
fundamental mode of the magnetic tube, T](t, z), namely 



i](t, z) = A(t)Wo(t, z) exp [iF(t, z)] exp 



Jo 



(W(T)dT 



(2) 



where A(f), Wo(t,z), and F{t,z) are real functions, and (L)(t) is 
the instantaneous frequency at time f. The function A(f) is pos- 
itive and represents the oscillation amplitude as a function of 
time. The functio n F(t,z) re presents a small phase shift that can 
be neglected (see iRudermaa201 lb) and, therefore, it is omitted 
hereafter The function Wo(t,z) satisfies the equation 



+ a/Wo = 0, 



(3) 



along with the boundary conditions Wo(t, +L/2) - due to 
photospheric line-tying. The function Wo(f, z) is defined so that 
Woit,z) > and max[Wo(f, z)] = 1. In Equation Q is the 
kink velocity, which in our model is defined as 



Ck(z, 



where 



Ckp, 
Cke, 



if z- < z < z+, 
otherwise, 



Ckp = 



2B2 



^ |io(pp +Pc) 



Cke = 



B2 



with fiQ the magnetic permeabihty of free space and 



Z± = ZO + UQt+ —, 



(4) 



(5) 



(6) 



with zo corresponding to the position of the center of the promi- 
nence thread with respect to the center of the magnetic tube at 
f = 0. Note that zo < if the thread is initially located on the left- 
hand side to the center of the tube, whereas zo > if the thread 
is initially located on the right-hand side. We also note that the 
flow velocity, mq, only explicitly appears in the expressions of z+ 
and Z-. 



The equatio n governing the a mplitude A(t) corresponds to 
Equation (66) of lRudermanI (1201 lb) , namely 



dt 
where 



nL/2 

(Pi+pJW^dz, 

J-L/2 



(7) 



(8) 



4 " 



L/2 



(Pi -pc)WoWn(r„)dz 



X 



J-L/2 



Pt(rn)-Pc 



Wow„(r„)dz. 



(9) 



-L/2 I An I 

In the expressions of / and E, pi is the internal density given by 



^ Pp, if z- < z < z+, 
' "c, otherwise. 



(10) 



and pt(r) is the density in the transitional layer, with pt(r„) the 
corresponding value at the n-th resonant position r - r„. The 
function w„ is the eigenfunction of the Alfven boundary value 
problem 



dz^ 



■ <^A,n('')'^n(r) = 0, 



(11) 



along with w„(r) = at z = ±L/2, where (r) = B^/jj.pt{r) 
is the local Alfven velocity squared. In Equation (fTTI) cjA,n{r) is 
the n-th Alfven eigenvalue. The eigenfunction w„ satisfies the 
normalization condition 



x: 



L/2 J 

-r-w„w„,dz = 6„ 

L/2C^ 



(12) 



where 6mn is the Kronecker delta. The position of the n-th reso- 
nance, r„, is given by the condition - w^jj(r„). Einally, the 
quantity A„ is given by 



A„ = - 



dr 



(13) 



We consider the case l/R - in o r der to compare with the 
expressions used bv lSoler & GoossensI (1201 ih . In the absence of 
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resonant damping F = and the equation for the amplitude A(t) 
(Equation ^) can be expanded to 



3.1. Finding Wo(t, z) and aj(t) 

To find the function Wo{t, z) we need to solve the boundary value 
problem defined in Equation Q with Wo(t, +L/2) = 0. In addi- 
tion, since the equilibrium is piecewise constant in z we need to 
provide additional boundary conditions at z - z±. Since z = z± 
correspond to contact discontinuities, the boundary conditions 
are th at Wait, z) is continuous at z - z± (iGoedbloed & Poedta 
12004 . Thus, the process here is equivalent to the process fol- 
lowed bv lSoler & Goos sens (201 j]) to find the i r func tion Q\(t,z). 
The reader is referred to .Soler & GoossensI (1201 ih for details. 
The expression for Wait, z) is 



Wo(f,z) = 



cos H- (^j if Z- < z < z+. 



(15) 



where ip - (p{t) is a time-dependent phase, and Di and D2 are 
given by 



D2 



(16) 



The expression for Wo(t,z) takes into account the condition 
max[Wo(f,z)] = 1. 

The time-dependent frequency, a>{t), requires to obtain first 
the time-dependent dispersion relation by providing additional 
boundary conditions which enable us to eliminate the time- 
dependent phase, (p. These additional conditions are that dWo/dz 
is continuous at z = z±. The expression for is 



Z- - arctan \ — - cot 

Ckp [ Ck< 



(f^cot[-(z_.^)l|, (17) 



and th e dispersion relation is Equation (14) of lSoler & GoossensI 
(1201 1^. The fundamental mode corresponds to the solution with 
the lowest frequency. By performing the first-order expansion of 
their Equation (14 ) with respect to the small parameter Pc/Pp, 
ISoler & GoossensI (|201 ih obtained an approximate expression 



for <jj(f) (their Equation (15)), which in the limit pdpp ^ 1 real- 
istic of prominences reduces to 



dA(t) 1 doj Id/ 

^^2^dr^^'>^27d7^« = «- ^''^ 

We compare Equatio n (fT4l) to Equation (20) of 
ISoler & GoossensI (1201 Ih . The onl y difference between the 
equation for the amplitude used by ISoler & GoossensI (1201 lb 
and the equation derived using the more general method by 
lRudermanr (l201 lb) resides in the term with t he tempo ral deriva- 
tiveof / in Equation (fT4l l. The equation of ISoler & GoossensI 
(I2OI ll) does not take into account the dependence of / on time. 
For the fundamental mode and slow flows, / may be roughly 
independent on time. This can b e checked once Wo(f, z) is 
found. In that case, the equation of ISoler & Goosseni (1201 ll) is 
an approximate description of the fun damental rnode ampli tude. 

Using the method developed by iRudermaiil (1201 1 alibi) and 
summarized in the previous paragraphs, the problem of study- 
ing damped standing oscillations of flowing threads is reduced 
to finding the functions Wo(f, z) and A(f), along with the expres- 
sion for the instantaneous frequency a>(t). This is done in the 
following subsections. 



2ckp ^JL/Lp 



yj[L -Lp){L+ |Lp) - 4 (zo + uotf 



(18) 



This is the expression for oj (t) used in the following computa- 
tions. In the absence of flow, mq = and co becomes indepen- 
dent of time and corresponds to the normal mode frequency (see 
ISoler etal.ll20r0l) . 

Since Wo(f,z) is known, we can compute the integral / 
(Equation (0). The result is 



2pcD] 



2p,Dl 



z- + Ljl 
2 



£!^sinfcz- + L/2) 

4w \ Cke 



+ (Pp +Pc) 

^2(1) 



2 

z+ 



Cke . l^U) 

— sm — (z+ 

4w \Cke 



■L/2) 



Ckp 



2a) 

sin I — z+ + 1 

Ckp 



Ckp 



-z- + 1 



(19) 



As before, we take advantage of the fact that the fundamental 
mode is the solution with the lowest frequency. We perform the 
first-order expansion of Equation fT% with respect to the small 
parameter Pc/Pp to get the approximate expression of / valid for 
the fundamental mode, namely 

I ~{pp+ pc) (z+ - Z-) = (pp + Pc) Lp. (20) 

Thus, we find that the integral / is time-independent in this ap- 
proxima tion. We can now check th e equation for the amplitude 
used by ISoler & GoossensI (1201 ll) in the absence of resonant 
damping, i.e., l/R - 0. Wh en / is c onstant in time. Equation ([T4l i 
reverts to Equation (20) of lSoler & Goossens (201 1). The refore , 
the equation for the amplitude used bv lSoler & GoossensI (|201 ll) 
is justified in this approximation. 

With the use of the approximate expression for / given in 
Equation (l20l i. we obtain the general equation for the amplitude 
when l/R + (Equation (|7|), namely 



(21) 



' (Pp+Pc)ip 

The next step is to compute F in order to integrate Equation (I2TI 1. 
3.2. Finding r 

Equation (|9]) for F takes into account the possibility of multiple 
resonances. From hereon we assume that there is only one reso- 
nance position at r = ri, and we set N = 1 in Eq uation (|9]l- This 
is the same assumption as in jRuderman' ('201 lb, Section 5.1). 
This assumption was also used by Soler etal. (2010) and was 
later checked to be correct by Arregui et al. (2011). In addition, 
we take into account that the density in the evacuated part of the 
tube, pe, is set equal to the coronal density, pc, so that l/R = in 
the evacuated region. Thus, we can change the limits of integra- 
tion in Equation (|9]i from [-L/2,L/2] to [z-,z+]. Also, since pp, 
Pc, and pt(ri) are constants, we can move them out of the inte- 
grals. The same is true for |Ai| since this quantity is independent 
of z. Finally, the expression for F becomes 



r - ^ (pp -pc)(pt(^i) -pc) ^^2 



2B^ R 



|Ai 



(22) 



4 



R. Soler et al.: Damped oscillations of flowing threads 



with 



A 



WoWi{ri)dz 



At the resonance position r - ri, so ap- - cix^ By com- 
paring Equations (O and (fTTI) we see that Wo and Wn(ri) are the 
eigenfunctions of the same eigenvalue problem when 



Pt(n) 



Pp +Pc 



(24) 



So, we can put w„(ri ) proportional to Wq, namely w„(ri ) - QWq, 
where Qis a constant of proportionality, and Equation ( l23b be- 



A 



z+ - Z- _^ 
2 2tij 



2fa; 

sm I — z+ + '■ 

Ckp 



2(1) 

Ckp 



-Z- + : 



(25) 



As before, we only take into account the fundamental mode 
and perform the first-order expansion of Equation ( l25T l with 
respect to Pc/Pp- We also use the normalization condition 
(Equation ( fT2l i) with m = n = 1 to find that Q ^ c\^^IL}J~ in 
this approximation for the fundamental mode. Thus, we obtain 
the approximate A as 



(26) 



Now we compute |Ai |. To do so we use again the resonant 
condition up- - aP^ j(ri) and Equation ( l24l i to write 



2 , , 2 4('-) 

Wai('-) = w —j— 



i Pp +Pi: 

' 2p,(r) 



(27) 



Hence the expression for |Ai| is 



\^^^oP 



Pp +Pc 

We next write 



5pt 



5r 



(28) 



dpx 



dr 



4 Z 



(29) 



where F is a factor that depends on the form of the transverse 
density profile. For example, F = for a linear profile and 
F = 2/7T for a sinusoidal profile with ri = 7?. The expression for 
I All becomes 



|Ai|=F 



uP Pp - Pc 
2 I Pp+Pc' 



(30) 



Finally, we substitute all these results in Equation ( l22l l and arrive 
at the expression for F, namely 



(31) 



3.3. Time-dependent amplitude A(t) 

We substitute Equation dSTT l in Equation ( 1211 1 to write the equa- 



(23) tion for the time-dependent amplitude A(t) as 



dA / 1 dw\ 

= - + A, 

df \^ 2w df / 



with 



1 I Pp-Pc 



2nF Rpp+pc 



(32) 



(33) 



Equation (l32T i can be easily integrated using the expression for 
CL)(t) given in Equation ( fTSl ). The result is 



A(0 - Ao 



(l - Lp) (l -I- iLp) - 4 (zo + uotf 



(L-Lp)(L+iLp)-4z2 



1/4 



(34) 



with the function ^(f) given by 




"(0 = —^h- 



arctan 



arctan 



2(zo + uot) 



^ ^J{L - Lp) [l + iLp) - 4 (zo + UQtf ^ 
2z() 



^(L-Lp)(L+iLp)-4z2 



(35) 



The first factor on the right-hand side of Equation ( l34l i. i.e., 
Ao, represents the initial amplitude at f = 0. The second fac- 
tor on the right-hand side of Equation (l34b . i.e., that with the 
square brackets, is the sa me as that found in Equation (22) of 
ISoler & GoossensI (1201 l b. This term accounts for the change of 
the amplitude due to the motion of the prominence thread along 
the magnetic tube. This contribution can cause either amplifica- 
tion or damping depending on the instantaneous position of the 
thread with respect to the center of the magnetic tube. This term 
is equal to unity in the static case, i.e., for mo = 0. The third 
factor on the right-hand side of Equation ( l34b is an exponential 
factor that accounts for damping due to resonant absorption. 

To recover the result in the static case without flow, we eval- 
uate the limit of Q(f) when mq 0, namely 



lim Q(f) = (Dot, 



(36) 



where wo the frequency given in Equation ( fTSb with mo - 0. 
In that case, the exponential factor on the right-hand side of 
Equation (l34l l becomes exp(-7Wof), so that we consistently 
revert to the static case studied by , Soler et al, (2010) and 
lArregui etan(l201lh . 

4. Parametric study 

Here we explore the impact of the various model parameters on 
the amplitude. In particular we focus on the effects of the thick- 
ness of the transitional layer, IjR, and the flow velocity, mo. We 
assume a linear variation for the density in the transitional layer 
Figure [2] shows A(r)/Ao as a function of f/Po, with Pq - 
271 1 cl)() the instantaneous period at f = 0. We plot the results for 
several values of IjR and zolL, and compare the amplitude ob- 
tained in the presence and in the absence of flow. We find that 
the amplitude in the presence of flow is substantially different 
from that for mo - when the transitional layer is very thin, so 
that damping due to resonant absorption is very weak. In this 
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Fig. 2. A(f)/Ao as a function of f/Po, with Pq the instantaneous 
period at f = 0. The numbers next to the various lines indicate 
the value of IjR used. The solid line is the complete result, while 
the dashed line is the equivalent result but with mq = 0. We have 
used (a) zqIL - -0.25 and (b) zo/L - 0.1. In all computations 
Lp/L = 0.1, Mo/cAp - 0.1, andpp/pc = 200. 

case the variation of the amplitude is mainly governed by the ef- 
fect of the flow, described by the factor with the square brackets 
on the right-hand side of Equation ( l34b . It is possible to obtain 
amplification of the amplitude at short times when zq/L < and 
the transitional layer is very thin (top lines in Fig.|2^). As l/R 
increases, the damping term, which is the exponential factor on 
the right-hand side of Equation ( l34b . becomes more important 
and the amplitude is almost the same in both cases mq = and 
Mo 0, i.e., the effect of the flow is then of minor importance. 
When zo/L > (Fig.|2j5) amplification is not possible regardless 
the value of l/R and damping is always observed. 

Figure [3] displays the thread displacement as a function of 
t/Po for different values of the flow velocity (the remaining pa- 
rameters are indicated in the caption of the Figure). We obtain 
a progressive phase shift between the solutions corresponding 
to different velocities. This is so because the instantaneous pe- 
riod, P{t) - 27r/w(f), is a function of the flow velocity accord- 
ing to Equation ( fTSl l. However, for the set of parameters used in 
Figure[3]the amplitude is weakly affected by the flow velocity. In 
particular, the value l/R = 0.1 used in Figure [3] is large enough 
for the amplitude to be dominated by resonant damping. 

5. Implications for prominence seismology 

The results presented in this pa per have implications f or seismol- 
ogy of solar prominences (see lArregui et alj|2012bh . Adopting 
resonant absorption as damping mechanism, the inversion tech- 
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Fig. 3. Prominence thread transverse displacement, r]{t)lri(), as 
a function of f/Po, with Pq the instantaneous period at f = 0, 
for various values of uo/cAp indicated within the figure. In all 
computations Lp/L - 0.1, l/R - 0.1, zo/L = -0.25, and Pp/pc = 
200. 

nique of iGoossens et al.l (l2008l) has been apphed to the case of 
prominence thread oscillations (see 'Arregui & BallesteJ 120101: 
ISoler et al. 2010). Due to the high density contrast of the promi- 
nence plasma with respect to the coronal plasma, a direct es- 
timation of l/R is possible when measures of both period and 
damping time are available. The value of l/R inferred by seis- 
mology may be used to test thread models, because the trans- 
verse inhomogeneity length scale is a crucial parameter for the 
energy balance of the prominen ce threads wi t h the surround- 
ing hot coronal plasma (see, e.g.. |Poio"gilll994t ICirigliano et al.l 
2004; Labrosse et al. 2010). The analytic equation used to infer 
l/R is (IGoossens et al.ll2008b 

/ P Pp+Pc P 

- = F— ^ ^ ^ F—, (37) 

R TDPp- Pc Td 

where P and td are the (constant) period and the damping time 
obtained after fitting to the observed oscillation a harmonic func- 
tion with amplitude proportional to exp (-f/rp). This expression 
does not take the flow of the prominence thread into account. 
In this Section we test the robustness of seismological estimates 
when flow is present and Equation dJTl l is used in combination 
with typical fitting methods for the oscillation amplitude. In ad- 
dition, the influence of background noise is also studied. 

5.1. Example 

We synthetically generate a signal aiming to represent a promi- 
nence thread transverse oscillation detected with a real instru- 
ment. We use the following set of parameters: Lp/L = 0.1, 
l/R = O.I, zo/L = -0.05, pp/pc = 200, and mq/cap = 0.1. 
A linear density profile in the transitional layer is used. The 
corresponding theoretical transverse displacement is displayed 
in Figure Ua) with a dashed line. The time series corresponds 
roughly to 6 oscillation periods. This is approximately the num- 
ber of periods observed in time series of real events (see, e.g., 
iNing et al.ll2009l) . To represent the limited cadence of the instru- 
ment we perform a temporal sampling of this signal using, ap- 
proximately, 33 points per period. For a typical period of 3 min, 
this corresponds to a cadence of about 5 s. This cadence is of 
the same orde r as that used in recent observations (see, e.g., 
iLin et al.ll200l . Later, we add to the sampled signal a randomly 
generated noise. In this way we try to account for the contri- 
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bution of the background to the total signal and/or to consider 
unspecified instrumental uncertainties. Noise is generated with 
the IDL function randomu. Here, we specify the amplitude of 
the background noise as percentage of the signal amplitude at 
t = 0, namely 



(%)„ 



X 100, 



(38) 



where //noise is the noise amplitude and 770 is the signal amplitude 
at f = 0. It is more frequent in observational astrophysics to use 
the signal-to-noise ratio, S IN, which is defined as the power ratio 
between the signal and the background noise. Assuming that the 
power is proportional to the square of the amplitude, the relation 
between (%)„oise and S /N is 



S/N : 



100 



(39) 



For example, (%)noise - 10% corresponds to S /N - 100. In 
the present example we use (%)noise = 50%, which is equivalent 
to S /N = 4. Finally, a smoothing algorithm is applied to the 
noisy data. The resulting signal is plotted in Figure with the 
solid line. Since the added noise is randomly generated, we can 
obtain different synthetic signals depending on the seed used in 
the random number generator. The synthetic data displayed in 
Figure EJ a) is the one used in this particular example. 

Figure ISb) shows a wavelet power spectrum 
(iTorrence & Compol Il998h of the synthetic signal displayed 
in Figure |4}a). For comparison we overplot the instantaneous 
period of the original data. We see that the wavelet spectrum 
recovers well the period of the original data. However, since 
the effect of the flow on the period is not very strong and the 
duration of the time series is limited to 6 periods, the variation 
of the period with time is not evident in the wavelet spectrum. 
The added noise also makes diflicult the analysis near the end 
of the time series when the amplitude is low. Based on this 
result, we believe that the effect of the flow on the period would 
probably go unnoticed for an observer using wavelet analysis. 

We turn again to the synthetic signal of Figure 13 a) and try 
a different approach. We fit to the synthetic data a typical expo- 
nentially damped harmonic function as 



/(f) - a\ cos (fl2f + ^3) exp i-a^f) , 



(40) 



where a\, 02, 03, and 04 are four parameters to fit. The fitted 
period, P, and damping time, td, are 



2n 

P = —, 

a2 



1 

Td = — ■ 

04 



(41) 



The fit is made using the IDL function curvef it. We perform a 
fit to the sampled data before and after adding random noise. The 
result of both fits is shown in Figure|4jc), where we also plot the 
original curve for comparison. First of all, we obtain a progres- 
sive phase shift between the different curves because the fitted 
function (Equation (l40l i) does not take into account the variation 
of the period with time. The constant period obtained from the 
fit to the noisy data is indicated by the horizontal dotted line in 
FigurelUb). This fitted period agrees well with the maximum in 
the wavelet spectrum, but overestimates the actual period at later 
times. 

More importantly, we also find different damping rates for 
the three curves. The damping of both fitted functions is stronger 
than the actual damping. In particular, the fitted curve to the 
noisy data has the strongest damping, i.e., the lowest amplitude. 




Original 

Fitted (no noise) 
Fitted (witti noise) 




3 
t/Po 



Fig. 4. (a) Original (dashed) and synthetic (solid) data used in 
the seismological test, (b) Wavelet power spectrum for the di- 
mensionless period, P/P^, corresponding to the synthetic signal 
displayed in panel (a). The white solid line is the original data 
instantaneous period, whereas the horizontal dotted line is the 
period obtained from the fitting method. The red solid line de- 
notes 99% of confidence level, (c) Comparison of the original 
(dashed) and fitted signals with noise (solid) and without noise 
(dotted). 



The difference in the damping of the original curve and that of 
the sampled signal without noise can be attributed to the effect 
of the flow, whereas in the sampled signal with noise we have 
to take into account the additional influence of noise. This has 
a direct consequence on the seismological estimation of l/R us- 
ing the fitted P and T£, in Equation ( |37] ). In the original data we 
used l/R = 0.1, while the inferred values are l/R = 0.13 for the 
sampled signal without noise and l/R - 0.16 for the sampled 
signal with noise. We have not computed the uncertainty of the 
fitted parameters and, therefore, we do not know the relative er- 
ror on the estimation of l/R. However, the results of this example 
suggest that the presence of both flow and noise adds significant 
uncertainty to the determination of l/R. 
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5.2. Statistical study 

Now, we perform a statistical study of the influence of flow and 
noise on the determination of l/R and its uncertainty. The pur- 
pose here is to know whether the results obtained in the example 
of the previous subsection are general or, on the contrary, they 
are particular to this example only. We use the same parameters 
as in the previous example but vary the value of the flow velocity 
and the percentage of maximum background noise. 

First we fix (%)„oise = 50% (S /N = 4). We follow the 
procedure explained before to generate a set of 10"* synthetic 
noisy signals for various values of uq. The various synthetic sig- 
nals are generated using different initial seeds in the random 
number generator As in the example, we fit an exponentially 
damped harmonic function (Equation ( l40t ) to each one of the 
synthetic signals and infer the value of l/R using Equation dJTl ). 
Subsequently, we perform a histogram with the full set of 10"* 
estimated l/R for each particular mq. By visual inspection, we 
determine that the histograms follow the Gaussian distribution. 
Therefore, we fit to the histograms a Gaussian function of the 
form 

g(l/R) = bo exp J- ^^^^1^ J , (42) 

with bo the height of the Gaussian, fi the mean value, and cr 
the standard deviation. We display in Figure |3 a) the histograms 
and their Gaussian fits corresponding to three values of the flow 
velocity. For mq = the Gaussian mean value matches the actual 
l/R. As the flow velocity increases, the mean value is shifted 
towards larger l/R, but the width of the Gaussian is roughly the 
same regardless of the flow velocity. 

On the other hand. Figure |5];b) shows equivalent histograms 
but for a fixed value of the flow velocity and three different 
(%)noise- We find that the width of the Gaussian is strongly af- 
fected by noise, i.e., the larger the background noise, the wider 
the Gaussian. The width of the Gaussian is related to the uncer- 
tainty of l/R. Thus, we consistently obtain that the uncertainty 
increases when noise is increased, as expected. However, the 
mean value is not affected by the percentage of noise. The effect 
of noise is explored in more detail in Figure |6l which shows the 
fitted yu (Figure |6^) and cr (Figure |6]3) as a function of the maxi- 
mum background noise percentage for three flow velocities. We 
find that /i is independent of noise and its shift with respect to 
the actual l/R is only determined by the flow velocity. For the 
set of parameters used in this test, the mean estimated l/R when 
flow is present is larger than the actual value. On the contrary, 
the effect of flow on cr is minor compared to the effect of noise. 
We obtain that cr is a approximately linear with the maximum 
noise percentage. 

Up to here we have obtained that the seismologically inferred 
l/R tends to overestimate the actual value. Our goal now is to 
determine whether this is the general tendency or whether this is 
affected by the choice of the model parameters. There are sev- 
eral parameters that may affect the estimation of l/R. In partic- 
ular, the role of zo, i.e., the position of the thread with respect 
to the center of the magnetic tube at the beginning of the oscil- 
lation, is worth being studied. To do so, we fix uo/c^p - 0.1, 
(%)n(,ise - 50%, and select a particular value of zo/L (the re- 
maining parameters are the same as before). For a given zo/L 
we repeat the previously explained procedure to generate the set 
of synthetic signals and the corresponding estimations of l/R. 
We fit the Gaussian function to the obtained histogram and com- 
pute fi and cr as functions of zo/L. These results are displayed in 
Figure]?] Regarding jj. (Figure]?^), we find that the estimated l/R 
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Fig. 5. (a) Normalized histogram of the seismologically esti- 
mated l/R for three values of the flow velocity (indicated within 
the panel) and a maximum background noise corresponding to 
50% of the initial amplitude, (b) Same as panel (a) but with 
Mo/cAp = 0.1 and three percentages of maximum background 
noise (indicated next to the lines). In both panels, the dashed 
lines are the Gaussian fits and the vertical dotted line is the ac- 
tual l/R. 



is larger than the actual value except when zo/l^ 5 -0.1. When 
Zo is negative, i.e., the thread is initially located to the left of the 
center of the magnetic tube, the effect of flow is to increase the 
amplitude until the thread reaches the center of the tube, mean- 
ing that flow opposes resonant damping. This increase of the 
amplitude compensates the overestimation of l/R by the fitting 
method, so that the mean estimated l/R is closer to the actual 
value. On the contrary, when zo/L > -0.1 the flow contributes 
to damping during most of the evolution. This causes an over- 
estimation of l/R. The behavior of cr (Figure|7j)) also shows an 
increasing trend with zo / L. 

Based on the results shown in this section, we conclude that 
flow and noise have two different effects on the estimation of l/R 
using typical fitting methods in combination with the theoretical 
inversion formula. An exponentially damped harmonic function 
has been used as a proxy to the actual oscillation. On the one 
hand, flow tends to shift the estimated l/R with respect to the 
actual value. The shift depends specifically on the flow velocity 
and on the position of the prominence thread at the beginning of 
the oscillation but, in general, an overestimation of l/R is found 
in the statistical analysis. On the other hand, as expected, noise 
adds uncertainty to the estimation. This noise-related uncertainty 
might be partially reduced if some method able to subtract the 
background noise is applied to the raw data. However, some un- 
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Fig. 6. (a) Mean value, and (b) standard deviation, cr, of the 
Gaussian distribution of seismologically estimated IjR as func- 
tions of the percentage of background noise. The different sym- 
bols correspond to different flow velocities (indicated within the 
panels). The horizontal dotted line in panel (a) shows the actual 
value of IjR. 

certainty intrinsically related to the effect of flow might remain 
even if the influence of noise is completely removed. 

Real data are frequently affected by gaps in the time series. 
This issue has not been taken into account in the present anal- 
ysis. The consideration of gaps would add more complexity to 
the analysis and, probably, would result in a new source of er- 
ror. Appropriate methods accoun ting for the effect of ga ps in the 
signals should be used (see, e.g.. lCarbonell et anil992h . This is 
worth being investigated in future works. 

6. Conclusions 

In this paper we have studied the joint effect of resonant absorp- 
tion and flow on the amplitude of standing kink waves in promi- 
nence threads. The present work extends the previous paper by 
ISoler & Goossens j ( 20 1 1 ) that did not take resonant damp ing into 
accoun t. We have followed the method developed by iRudermanI 
(l20TTallbh to obtain a general analytic expression for the kink 
mode amplitude as a function of time which includes the effects 
of both resonant absorption and flow. 

We find that flow and resonant absorption can either be com- 
peting effects on the amplitude or can both contribute to the 
damping of the kink mode depending on the instantaneous po- 
sition of the dense thread within the prominence magnetic tube. 
For fast flows and short transverse inhomogeneity length scales 
the amplitude profile deviates from the classic exponential func- 
tion for resonantly damped kink modes in static flux tubes. From 



Fig. 7. (a) Mean value, and (b) standard deviation, cr, of 
the Gaussian distribution of seismologically estimated IjR as 
functions of zqIL. A maximum background noise of 50% and 
uo/cAp = 0.1 are used. The horizontal dotted line in panel (a) 
shows the actual value of l/R. 



the observational point of view, to determine the location of the 
dense plasma within the magnetic tube might be difficult since 
the footpoints of the magnetic tube are not seen in the observa- 
tions. 

The implications of our results for seismology of solar 
prominences have been explored. We have test the robustness of 
seismological estimates of the transverse inhomogeneity length 
scale. We have used synthetic data aiming to mimic real obser- 
vations and have performed a statistical study. Our results show 
that the presence of flow can significantly affect the estimation 
of the transverse inhomogeneity length scale. Statistically, we 
find that this parameter is overstimated when an exponentially 
damped harmonic function, which does not take flow into ac- 
count, is used to fit the actual oscillation. The presence of ran- 
dom background noise and/or intrumental errors adds further un- 
certainty to this estimation. 

The seismology of flowing prominence threads using 
damped kink waves is more challenging than that of their static 
counterparts because of the effect of flow on the amplitude. 
The presence of flow adds complexity to the behavior of the 
oscillations and has a direct impact on seismology (see, e.g., 
Terradas et al. 201 1). Caution needs to be paid to the seismolog- 
ical estimates that do not take the influence of flow into account. 
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